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ABSTRACT 

Aims. This paper presents a detailed comparison between high-redshift observations from the VIMOS-VLT Deep Survey (VVDS) and predictions 
from the Munich semi-analytical model of galaxy formation. In particular, we focus this analysis on the magnitude, redshift, and colour distribu- 
tions of galaxies, as well as their clustering properties. 

Methods. We constructed 100 quasi-independent mock catalogues, using the output of the semi-analytical model presented in De Lucia & Blaizot 
(2007). We then applied the same observational selection function of the VVDS-Deep survey, so as to carry out a fair comparison between models 
and observations. 

Results. We find that the semi-analytical model reproduces well the magnitude counts in the optical bands. It tends, however, to overpredict the 
abundance of faint red galaxies, in particular in the V and z' bands. Model galaxies exhibit a colour bimodality that is only in qualitative agreement 
with the data. In particular, we find that the model tends to overpredict the number of red galaxies at low redshift and of blue galaxies at all 
redshifts probed by VVDS-Deep observations, although a large fraction of the bluest observed galaxies is absent from the model. In addition, the 
model overpredicts by about 14 per cent the number of galaxies observed at 0.2 < z < 1 with I AB < 24. When comparing the galaxy clustering 
properties, we find that model galaxies are more strongly clustered than observed ones at all redshift from z = 0.2 to z = 2, with the difference 
being less significant above z — 1. When splitting the samples into red and blue galaxies, we find that the observed clustering of blue galaxies is 
well reproduced by the model, while red model galaxies are much more clustered than observed ones, being principally responsible for the strong 
global clustering found in the model. 

Conclusions. Our results show that the discrepancies between Munich semi-analytical model predictions and VVDS-Deep observations, particu- 
larly in the galaxy colour distribution and clustering, can be explained to a large extend by an overabundance of satellite galaxies, mostly located 
in the red peak of the colour bimodality predicted by the model. 

Key words. Cosmology: observations - Cosmology: large scale structure of Universe - Galaxies: evolution - Galaxies: high-redshift - Galaxies: 
statistics 



1. Introduction 

Most accurate cosmological probes tend to favour a flat KCDM 
cosmological model and strengthen the hierarchical growth of 



structure scenario (e.g. Riess et al. 


1998; Spergel et al. 2003 


Tegmark et al. 2004 Komatsu et al. 


2009). In this picture, it is 



believed that the onset of galaxy formation arises inside dark 
matter haloes and that the cosmic history of galaxies follows the 
hierarchical evolution of haloes. 
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Despite the undeniable successes of the ACDM model to 
explain a broad variety of astrophysical observations, the de- 
scription of stellar mass assembly and star formation activity 
in this framework and its confrontation to observations remain 
challenging (e.g. De Lucia 2 009 |l. The most recent deep galaxy 



spectroscopic surveys (e.g. 



2003 Le Fevre et al. 2005a 



Davis et al~ | 2003| Dickinson et al 



Lilly et al.|2007 1, that enabled the 



census of galaxy properties up to z — 2, have confirmed that the 
so-called downsizing scenario ( |Gavazzi et al.|1996{|Cowie et al.| 
1996), is a primary feature of galaxy formation. About half of 
the stellar mass in massive galaxies observed at the present time 
was already in place at z — 1 (e.g. Bundy et"al~||2005 Cimatti 



[eTaTlpOOfJl |Arnouts et a!J|2"007l |Ilbert et al.||2010| >, while the 
number density of less massive galaxies continues to rise with 
cosmic time even for redshifts below unity. Most massive galax- 
ies cease to efficiently produce new stars earlier than lower mass 
galaxies, which keep producing stars until very recent times (e.g. 
Tres se et al.|2007| l. This scenario is also supported by most re- 
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cent observations o f the VIMOS-VLT Deep Survey (VVDS,|Le| 
Fevre et al.|2005a) , one of the largest spectroscopic surveys of 



tering strength than observed, suggesting a significant deficit of 



distant galaxies. In particular, the VVDS-Deep sample provides 
one with a description of the high-redshift Universe from z — 0.2 
to z — 2, where the evolution of the galaxy luminosity function 
( Ilbert et al.|20 05 ; Zucc a et al.|2006| l, galaxy stellar mass func- 
tion (jP ozzetti et al. 2007), star formation rat e (|Tresse et al. 2007 ; 
Bard elli et al.||2009| l, " colour bimodality (|Franzetti et al.||2007) , 
and galaxy clustering ( Le Fevre et al.|2005b) cannot be only ex- 
plained by a simple scenario of hierarchical growth of baryons. 

The complex physical processes involving baryons on galac- 
tic scales likely explain why the physics of stellar mass as- 
sembly and star formation is still puzzling in the KCDM 
model. Fortunately, large N-body simulations coupled with 
semi-analytical treatments of baryonic processes or fully hydro- 
dynamical simulations have become sufficiently sophisticated to 
draw detailed and reliable predictions for galaxy properties (e.g. 
IWeinberg et aTpXMl |Springel et aLpXB) |Kim et al.| |2009b). 
Current semi-analytical models can be tuned to optimally re- 
produce some of the basic galaxy properties observed in the 
local Universe such as the galaxy luminosity function or the 
Tully-Fisher relation. Galaxy mock samples constructed from 
these models are widely used for a large variety of purposes. 
However, it is not yet demonstrated that numerical simulations 
together with semi-analytical model can successfully describe 
galaxy properties at high redshift and the physics driving the 
evolution of the various galaxy populations. There have been few 
detailed comparisons made between their predicted galaxy prop- 
erties and a broad range of observational measurements at high 
redshift. Most comparisons were focused on single observations 
that has provided important clues, but still a limited perspective 
on the global galaxy evolution scenario. 



Bower et al. ( 2006 1 compare the __-band luminosity function, 



galaxy stellar mass function, and cosmic star formation rate from 
the Durham model (Bow er et al.|20 06 ) with high-redshift obser- 
vations. They find that their model match the observed mass and 



luminosity functions reasonably well up to z — 1 . Kitzbichler & 
| White| ( p007| > compare the magnitude counts in B,R, I, K bands, 
redshift distributions for /_-band selected samples, B- and K- 
band luminosity functions, and galaxy stellar mass function from 
the Munich model ( |Croton et al.|2006||De~Lucia & Blaizot|2007| l 
with deep surveys measurements. They find that the agreement 
of the Munich model with high-redshift observations is slightly 
worse than that found for the Durham model. In particular, they 
find that the Munich model tends to systematically overestimate 
the abundance of relatively massive galaxies at high redshift. 
They note, however, a non-negligible dispersion between vari- 
ous observations. 

The predicted galaxy clustering properties by semi- 
analytical models have been mainly compared with local ob- 
servations (e .g. |Springel et aL|2"0"0"5j |Li et al.||2006] |2007| |Krm 



et al.|2009a] ). Similar conclusions are obtained when confronting 
mode ls to 2dFGRS ( |Colless et al.|2001| > and SDSS ( [York et al.| 
2000) measurements: semi-analytical models are able to repro- 



duce the overall observed clustering properties, but fail to pre- 
dict the measurements in details, e.g. the different clustering 



of red and blue galaxies (Springel et al. 2005). For the high- 
redshift Universe, only few data-model comparisons have been 
carried out ([McCracken et al.||2007) |Coil et al.pOOSI |Meneux| 



et al.|2008 2009 1. These studies highlight a number of discrep- 
ancies between model predictions and observations. In partic- 



ular, Coil et al. (2008]) show that red model galaxies are more 
strongly clustered than observed at z — 1, particularly on small 
scales. Blue galaxies in the model show instead a lower clus- 



blue satellite galaxies in the semi-analytical model of Croton 



et al. (20061. They interpret these discrepancies as due to an 



incorrect modelling of the colours of satellite galaxies at high 
redshift. 

In this paper, we compare the basic galaxy properties ob- 
served in the VVDS-Deep spectroscopic sample at 0.2 < z < 2 
to the outputs of the Munich model. This is the first of a series of 
papers that aim at a detailed comparison between model predic- 
tions and a complete set of galaxy properties and their evolution 
with cosmic time, derived from the VVDS-Deep observations. 
In this first paper, we focus on the magnitude counts, redshift 
distribution, colour bimodality, B-band and /-band luminosity 
distributions, and the global clustering measurements. The de- 
tailed dependence of galaxy clustering on luminosity and intrin- 
sic colour is the subject of a forthcoming paper (Me neux et al.| 
|2010| hereafter paper II). 

In Sect. 2 we summarise the basic characteristics of the 
VVDS-Deep sample and the Munich semi-analytical model. In 
Sect. 3 we compare the observed magnitude counts, redshift 
distribution, colour bimodality, B-band, and /-band luminosity 
distributions with model predictions. In Sect. 4 we present the 
galaxy clustering comparison and we discuss our results in Sect. 
5. Throughout this paper we assume a flat AC DM cosmology 
with £_ M = 0.25, __ A = 0.75 and H = 100 h km ■ s 1 ■ Mpc" 1 . 
All magnitudes are quoted in the AB system and for simplicity 
we denote the absolute magnitude M x — 5 log (h) as M x . 



2. VVDS-Deep sample and Munich galaxy formation 
model 



2.1. The VVDS-Deep spectroscopic sample 



The VIMOS-VLT Deep Survey (VVDS, |LeFevre et al.|20 05a) 
is a deep spectroscopic survey of galaxies performed with the 



Visible Multi-Object Spectrograph (VIMOS, |Le Fevre et al. 



2003) at the European Southern Observatory's Very Large 
Telescope (ESO-VLT). The deep part of the survey (VVDS- 
Deep) spans the apparent magnitude range 17.5 < Iab < 24 and 
contains about 8, 700 galaxies over an area of 0.49 deg 2 . In ad- 
dition to VIMOS spectroscopy, the field has a multi-wavelength 
photometric coverage: B, V, R, I from the VIRMOS Deep 
Imaging Survey (VDIS, |Le Fevre et al.|2004| ), u*, g', r', i' , z' 
from the Canada France Hawaii Telescope Legacy SurvejQ 



(CFHTLS-D1, Goranova et al. 



20091 [Coupon et al. 200 9), and 



a partial coverage in J and K bands from VDIS ( Iovino et al.| 
|2U03| ). 

Unless otherwise stated, in this analysis we use only the 
galaxies with secure redshift, i.e. which have a redshift confi- 



dence level greater than 80% (flag 2 to 9, see Le Fevre et al 



2005a i. This provides us with a spectroscopic sample of 6,582 
galaxies. The accuracy of the redshift measurements is of cr z = 
9.2 ■ 10~ 4 . The observational spectroscopic strategy allows us 
to reach an average spectroscopic sampling rate of 27% across 
the field. The observations, survey strategy, and basic proper- 
ties of the VVDS-Deep sample are described in detail in |Le| 
Fevre et al. (2005a). The VVDS-Deep spectroscopic catalogue is 
publicly available through the CENter for COSmology database 
(CENCOS) sit<5 

Absolute magnitudes for VVDS-Deep galaxies have been 
computed by Ilbert et al. ( 2005 1 and properly k-corrected us- 



http://terapix.iap.fr/cplt/table_syn_T0006.html 
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ing the spectral template that best fits the B, V, R, I photometry. 
The VVDS-Deep apparent and absolute magnitudes used in this 
study are not corrected for intrinsic dust attenuation. 

2.2. The galaxy formation model 

In this study, we take advantage of the semi-analytical model 
presented in De Lucia & Blaizot ( 2007| l, whose outputs are 
publicly available through the Millennium Simulation database 
site^j The model is implemented on the Millennium Run N- 
body simulation, a large dark matter N-body simulation which 
traces the hierarchical evolution of 2160 3 particles between 
Z = 127 and z = in a cubic volume of 500 3 /i~ 3 Mpc^ . 
It assumes a KCDM cosmological model with (O m = Qrf m + 
Cl b , Q h , Q/„ h, n, erg) = (0.25, 0.75, 0.045, 0.73, 1, 0.9). 
The resolution of the N-body simulation, 8.6 x 10 8 hT l M Q , cou- 
pled with the semi-analytical model allows one to resolve haloes 
containing galaxies with a luminosity of 0.1L* with a mini- 
mum of 100 particles. Haloes and sub-haloes are identified from 
the spatial distribution of dark matter particles using a stan- 
dard friends-of- friends algorithm and the SUBFIND algorithm 
( Springel et al.|200l) . All sub-haloes are then linked together to 
construct the halo merging trees which represent the basic in- 
put of the semi-analytical model. Details about the simulation 



and the merger tree construction can be found in Springel et al. 
( |20Q5] >. 

Galaxies are simulated on top of the dark matter simula- 
tion using the information contained in the halo merging trees. 
The semi-analytical model (SAM) used in this study includes 



ingredients and methodologies originally introduced by White 
& Frenk| ( 1991 1 and later refined by Kauffmann & Haehnelt 
(|2000|i, ISpringel et al.| (pOOT), |De Lucia et al.| ( |2004| i, |Croton 



|et al.| ( |2606) , and |De Lucia & Blaizot| ( |2007| >. The model includes 
prescriptions for gas accretion and cooling, star formation, feed- 
back, galaxy mergers, formation of super-massive black holes, 
and treatment of "radio mode" feedback from galaxies located 



at the centres of groups or clusters of galaxies (see Croton et al 
2006 for details). One important aspect in the modelling of 
galaxy properties at high redshift is dust extinction. This ef- 
fect is particularly critical when modelling galaxy magnitudes 
( [Kitzbichler & White|2007HFontanot et al.|2009b|). In this work 
we have used the dust extinction model of IDe Lucia & Blaizotl 
(2007), where dust extinction is parametrised using a homoge- 
neous interstellar medium component combined with a simple 
model for molecular clouds attenuation around newly formed 
stars. In the following, we will often refer to this model as the 
"Munich" model. 



2.3. Mock samples construction 

In order to compare the galaxy properties predicted by the model 
with VVDS-Deep observations, we generated mock samples that 
cover the same volume probed by the VVDS-Deep sample, as 
viewed by an observer at z = 0. These have been constructed us- 
ing the MoMaF facility (Blai zot et al.|2 005), which converts the 
output of a galaxy formation model into simulated catalogues of 
observations. In particular, the light-cone construction adopted 
in MoMaF avoids replication and finite-volume effects in con- 
structing mock survey catalogues. 

We built 100 quasi-independent mock samples of 0.49 deg 2 
providing the apparent magnitudes in the VDIS and CFHTLS 
filters for all galaxies. In computing magnitudes, we use the ap- 



http://www.mpa-garching.mpg.de/millennium/ 



propriate filter transmission curves of the VDIS and CFHTLS 
photometric bands. These mock samples are linked to the main 
SAM output, providing the full set of intrinsic galaxy properties 
that the model calculates, as for instance galaxy absolute mag- 
nitudes and host halo properties. The VVDS-Deep galaxy selec- 
tion criterion, 17.5 < Iab < 24, has been applied to the mock 
samples. These mock samples mimic the expected VVDS-Deep 
sample with a spectroscopic sampling rate of 100% and a uni- 
form sampling on the sky. Hereafter, we will refer to them as the 
"complete" mock samples: Cmocks. 

From the Cmocks we constructed a second set of mock sam- 
ples that include the detailed observational selection function 
and biases of the VVDS-Deep sample, and in particular, its com- 
plex angular sampling. These have been constructed using the 
following procedure: 

Photometric mask. In the VDIS images, objects falling in bad 
photometric regions or in areas where the presence of bright 
stars saturates the CCD, have been removed from the pho- 
tometric catalogue. Consequently, no objects inside these re- 
gions have been targeted for spectroscopy. In practice, these 
regions have been identified and coded into a photometric 
mask, that we applied to the mock samples. This repro- 
duces in the simulated catalogues the various empty regions 
present in the angular distribution of objects in the VVDS- 
Deep (see Fig. [2J. 

Target sampling rate. The SSPOC software (Bot tini et al.| 
|2005| l has been used to prepare VVDS-Deep observations, 
performing the design of the slit masks for the planned 
VIMOS pointings. This software, which accounts for the 
precise shape of the VIMOS field-of-view (4 quadrants de- 
limited by an empty cross), enables one to optimise the 
slit positioning and to maximise the number of objects ob- 
served in spectroscopy from a given list of potential targets. 
Similarly, we applied SSPOC to the mock samples giving as 
input the list of VIMOS observed pointings in the VVDS- 
Deep. SSPOC requires the angular size of objects. Since this 
property is not attributed to mock galaxies by the model, we 
randomly assigned an apparent radius to mock galaxies, in a 
way to reproduce the observed distribution of apparent radii 
as a function of the selection magnitude in the VVDS-Deep. 
In the optimisation process, SSPOC tends to preferentially 
select objects with smaller apparent angular size. This intro- 
duces a mild bias against large galaxies, on average brighter, 
which we will refer to as the target sampling rate (TSR) in 
the following ( |Ilbert et al.|2005| >. 

Spectroscopic success rate. Galaxy redshifts in the VVDS- 
Deep sample have been determined from the observed spec- 
tra. In this process, a confidence class is given to each tar- 
geted object to quantify its type and the confidence level 
on the redshift determination ( |Le Fevre et al.||2005a| l. In 
general, redshifts of galaxies with brighter apparent magni- 
tudes, and in turn with higher signal-to-noise spectra, have a 
larger probability to be measured. The success in determin- 
ing galaxy redshifts is thus a function of the apparent mag- 
nitude. This introduces a bias against faint objects which we 
will refer to as the spectroscopic sampling rate (SSR) in the 
following (Ilbert et al. 2005). To reproduce this selection bias 
in the mock samples, we randomly assigned a flag to each 
mock galaxy, in a way to match the dependence of the frac- 
tion of the different redshift confidence flags on the apparent 
magnitude in the VVDS-Deep. We kept in the mock samples 
the galaxies with flag 2 to 9, corresponding to objects with 
secure redshift in the VVDS-Deep (confidence level greater 
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Table 1. Definition and properties of VVDS-Deep samples. 
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Fig. 1. Fraction of the number of galaxy with secure redshift 
(flag 2 to 9) over the total number of spectroscopic targets in 
the VVDS-Deep as a function of the / apparent magnitude (his- 
togram). The solid curve is the average fraction over the Omocks 
samples. 



than 80%, Le Fevre et al.|2005a l. Fig. [T] shows the fraction 
of galaxies with secure redshift measurements over the to- 
tal number of spectroscopic targets in the VVDS-Deep, as 
a function of the apparent magnitude of selection. The solid 
curve in the same figure shows the mean measurement from 
the mock samples. 

By applying this procedure to the Cmocks samples, we end 
up with a new set of mock samples having on average a mean 
sampling rate of 27% (identical to that of the VVDS-Deep), and 
that include the detailed selection function and observational bi- 
ases of the VVDS-Deep spectroscopic sample. The angular dis- 
tribution of galaxies in the mock sample 1 obtained following 
this procedure is shown as an example in Fig. [2] and compared 
to that of the VVDS-Deep secure redshift sample. In the follow- 
ing, we will refer to these mock samples as the "observed" mock 
samples: Omocks. 

3. Number counts 

3.1. Magnitude counts 

We compare the galaxy magnitude counts in the VDIS /-band 
(selection band) and in the CFHTLS u* , g' , r' ,i' ,z' bands with 
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Fig. 2. Angular distribution of galaxies in the VVDS-Deep (left) 
and in the Omocks sample 1 (right). The galaxies are shown with 
the dots and the background coloured areas encode the surface 
density of objects. The solid contours encompass the regions dis- 
carded by the photometric mask. 



those predicted by the SAM. In this comparison, we use all 
galaxies with VDIS photometry for /-band counts, while for 
u* ,g' , r', i',z' counts, we use all galaxies with CFHTLS-D1 pho- 
tometry falling in the VDIS field. The effective areas probed are 
1.19 deg 2 and 0.89 deg 2 respectively. We used Cmocks sam- 
ples of same angular sizes. The VDIS and CFHTLS photomet- 
ric catalogues reach limiting magnitudes (80% completeness for 
point-like sources) of 24.6, 26.2, 25.9, 25.4, 25.1, 24.6 respec 
tively for /, u*,g', r',i', z' bands (see |McCracken et"aLl|2003 



Goranova et al. 2009). When computing the magnitude counts in 



the VVDS data, we carefully excluded stars using colour-colour 
diagrams following McCracken et al. (2003 1. 

We first compare the /-band galaxy counts in the VDIS field 
with published measurements from large and deep photometric 
surveys. This comparison is shown in Fig. [3] We find a very 
good agreement between the VVDS and other published /-band 
counts in the literature (see also McCracke n et al.|2003) >. When 
comparing to the SAM predictions, we find that the predicted 
/-band counts match very well the VVDS-Deep ones at / > 20. 
The model, however, slightly underpredicts the number of galax- 
ies brighter than / ^ 20. On average, the mock samples con- 
tain 10.3% less galaxies than observed at 17.5 < / < 24 in 
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the VVDS-Deep. In Fig.|4]we compare the predicted magnitude 
counts in the CFHTLS bands with the observational measure- 
ments. Note that in this figure, the CFHTLS counts extend above 
24 mag, while the SAM mock samples have been cut at / = 24 
to match the VVDS-Deep selection function. We find that the 
observed and predicted counts are in very good agreement, ex- 
cept in the i'- and z'-bands where the SAM predicts more galax- 
ies at the faintest magnitudes. While in the /'-band the counts 
are compatible within the 3<x field-to-field dispersion among the 
mock samples, the z'-band counts deviate from the VVDS-Deep 
at more than 3<x for magnitudes fainter than z' = 22. 

Similar small discrepancies between model predictions and 
observed magnitude counts were found by Kitzbichler & White 
(2007) who, however, used observational measurements based 
on smaller samples. In particular, they compared the counts in 
the B, R, I bands as observed over an area of 0.2 deg 2 in the 
Hubble Deep Field North (HDF-N), and in K-band, in fields 
of area smaller than 0.1 deg 2 . They found a very good agree- 
ment between model predictions and observations in the B, R, I 
bands, with only a small underestimation of the predicted counts 
at / < 19.5. For ^T-band counts, the agreement they found was 
less good with the counts becoming discrepant at K > 21 by a 
factor of 2. These results are fully consistent with ours. 

3.2. Redshift distribution 

To provide a fair comparison of the VVDS-Deep redshift distri- 
bution with that predicted by the SAM, we first consider only the 
VVDS-Deep galaxies with secure redshifts (flag 2 to 9), which 
represent 75% of the full galaxy sample. The remaining 25% of 
the sample is made of galaxies with a redshift confidence level 
of about 50% (flag 1, 17%) and unclassified objects for which no 
redshift measurement has been possible (flag 0, 8%). Note that 
this latter class possibly includes stars. As redshifts are available 
only for about 27% of the galaxies with 17.5 < I < 24, we nor- 
malise the redshift distribution to the total number of galaxies 
in the parent photometric catalogue. The corresponding redshift 
distribution is shown in Fig.[5]with the solid thin histogram. This 
normalisation assumes that the spectroscopic sample is statisti- 
cally representative of the parent photometric catalogue and that 
there are no biases introduced by the particular spectroscopic 
targeting strategy adopted. In reality, both these assumptions are 
incorrect. 
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Fig. 3. Comparison of the magnitude counts in the VDIS /-band 
with SAM predictions and previously reported measurements 
by|Postman et al.| ( 1998) > in the KPNO survey, Leauthaud et al. 
( [20071 1 in the COSMOS- ACS field, and |Metcalfe et al.|p001] > in 
the Hubble Deep Field-North & South. The VVDS-Deep counts 
are shown with the empty squares and the error bars correspond 
to their associated Poisson errors. The solid curve and associated 
shaded areas are respectively the mean and the l<x, 2<x, and 3cr 
field-to-field dispersions among Cmocks samples. 



In order to account for the observational biases and estimate 
their influence on the observed redshift distribution, we use the 
information of the survey target sampling rate (TSR) and spec- 
troscopic sampling rate (SSR) available for each galaxy of the 
sample. These two functions have been originally defined to 
correct for the observational biases in the measurement of the 



galaxy luminosity function ( Ilbert et al.|2005|[Zucca et aL 2006 ). 
The SSR accounts for the fraction of objects without a reliable 
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Fig. 4. Comparison of the magnitude counts in the CFHTLS u* , g', r', i', z' optical bands with SAM predictions for 17.5 < I <24 
galaxies. The CFHTLS counts are shown with the squares and the error bars correspond to their Poisson errors. The curves and 
associated shaded areas are respectively the mean and the lcr, 2<x, and 3<x field-to-field dispersions among the Cmocks samples. 
SAM predictions are for samples explicitly cut at / = 24, while CFHTLS counts extend above this limit. All galaxies with / < 17.5 
have been removed from the counts. 



redshift determination (flag and flag 1 galaxies) and corrects 
for the fact that the success rate of redshift measurements de- 
creases at fainter apparent magnitudes. The TSR instead, cor- 
rects for the selection biases introduced by SSPOC software in 
the spectroscopic mask preparation (Bottini et al. 2005| l, e.g. 
its tendency to target objects with small apparent angular size. 



Within the Ultra-Deep part of the VVDS survey (Le Fevre et al. 



2010), a small fraction of flag 0, flag 1, and flag 2 galaxies have 
been reobserved. These reobservations, representing 4% of the 
total number of objects in the spectroscopic catalogue, have per- 
mitted us to refine the measurement of the TSR and SSR func- 
tions, allowing us to better statistically account for the spectro- 
scopic incompleteness. The VVDS-Ultra-Deep observations, as 
well as details on the calculation of the TSR and SSR are given 
in |Le Fevre eta^pOTol ). 

We correct the raw redshift counts, including flag 1 galax- 
ies, by weighting each galaxy by w = (TSR x SSR) -1 . The 
connected redshift distribution is shown in Fig. [5] with the solid 
thick histogram. The solid curve and associated shaded area cor- 
respond to the predicted mean N{z) and the lcr field-to-field dis- 
persion among Omocks samples, which provides an estimate of 
the sample variance in the model. The predicted mean N(z) of 
the Cmocks samples is plotted with the dotted curve. 

We first note that VVDS-Deep observational biases have lit- 
tle effect on the shape of the predicted redshift distribution in 
the mock samples: the mean N(z) obtained from Omocks is very 



similar to that of Cmocks. The SAM well reproduces the shape 
of the observed distribution between z — 1 and z — 1 .8 while it 
shows an excess of about 14 per cent with respect to the observa- 
tions at 0.2 < z < 1. Similar trends were found in Kitzbic hler &| 



White (2007 1 who noted that the predicted redshift distribution 
was higher than observed over the redshift range 0.5 < z < 1 .5 
for K < 21.8 galaxies. As we will discuss in the next section, 
this excess is related to an excess of red galaxies at z < 1 in the 
model. Furthermore, we find that the model does not account for 
the tail of the distribution at redshifts above z — 2 and that ex- 
tends to z - 4 in the VVDS-Deep sample. In fact the SAM does 
not predict any galaxy at z > 3 with / < 24. The dip observed 
at 1.8 < z < 3 in the VVDS-Deep uncorrected N(z) is a purely 
observational effect usually referred as the "redshift desert". It 
is due to the lack of spectral features in galaxy spectra observed 
within the wavelength window function of the spectrograph at 
these redshifts ( |Le Fevre et~aT1|2005a| l. VVDS-Ultra-Deep re- 
observations, based on spectra measured on a larger wavelength 
range (VIMOS LR-Blue plus LR-Red grisms), allow us to cor- 
rect for this observational effect and to repopulate the "redshift 
desert" at 1.8 < z < 3. Indeed, a large fraction of the reobserved 
galaxies falls in this part of the redshift distribution, and thus by 
properly weighting the original flag 1 to 9 galaxies we are able 
to statistically account for the missing fraction of objects at these 
redshifts (see Le Fevre et al., 2010). 
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Fig. 5. Comparison of the VVDS-Deep with the SAM redshift 
distributions for 17.5 < Iab < 24 galaxies. The VVDS-Deep 
N(z), while including only flag 2 to 9 galaxies, is shown with 
the solid thin histogram. The thick histogram corresponds to the 
corrected N(z) in the VVDS-Deep when accounting for the ob- 
servational biases of the survey (see text for details). The solid 
curve and associated shaded area are the mean and the l<x field- 
to-field dispersion among the Omocks, while the dotted curve is 
the mean prediction of the Cmocks. 



3.3. Galaxy luminosity and intrinsic colour distributions 

Galaxies are found to have a bimodal colour distribution and 
rest-frame colours are commonly used to differentiate red mas- 
sive early-type from blue star-forming late-type galaxies. Here 
we define these two populations on the basis of the rest-frame 
B - I colour distribution. We made this particular choice be- 
cause it corresponds to the optimal rest-frame colour that can be 
measured both in the VVDS-Deep and in the model (the rest- 
frame {/-band is not available in the SAM). We considered six 
redshift intervals from z - 0.2 to z = 2.1, and compare the rest- 
frame fi-band, /-band, and B - I distributions from the model 
with the observational measurements. The definition of the sam- 
ples and their basic properties are given in Table [T] and Table [2] 
The galaxy luminosity and rest-frame colour distributions (nor- 
malised to unity) in the VVDS-Deep and the SAM are presented 
in Fig. [6] and Fig. [7] 

The predicted rest-frame B-band distributions in the differ- 
ent redshift intervals are in good agreement with VVDS-Deep 
observations (the apparent discrepancy in the peak height in the 
last high-redshift interval is not statistically significant given the 
small number of galaxies involved). This indicates that the B- 
band luminosity function may be well reproduced in the SAM 
up to z — 2 as also found by Kitzbichler & White ( |2007| ) at 
0.2 < z < 1.2. Instead, we find that the rest-frame /-band dis- 
tributions are systematically skewed towards bright magnitudes 
in the SAM. The effect is particularly important at z > 0.9 
where the model tends to significantly overestimate the number 
of galaxies with -23 < Mj < —21. A similar trend was found in 
Kitzbichler & White ( 2007 ) when confronting the /if -band lumi- 
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Fig. 6. Comparison of the rest-frame fi-band (right panels) and I- 
band (left panels) distributions in six magnitude-limited samples 
from z - 0.2 to z - 2.1. The histograms are VVDS-Deep mea- 
surements while the curves and associated shaded areas corre- 
spond to the mean and the lcr field-to-field dispersion among the 
Omocks samples. The dashed curves are the mean distribution 
in the Omocks without including dust extinction in the model, 
while the dotted ones are the mean predictions of the Cmocks. 



nosity function, although the small observational samples they 
used did not allow them to reach firm conclusions. 

When comparing the intrinsic colour distributions, we find 
that the rest-frame B — I colour distribution is clearly bimodal 
both in the VVDS-Deep and in the model, but the agreement is 
only qualitative. As shown in Fig. [7] the SAM does not repro- 
duce quantitatively the observed intrinsic colour distributions: 
there are many fewer very blue galaxies (i.e. with B - I 0.3) 
and much more "green valley" galaxies (i.e. with B - I 1) 
in the model than in the observations, at all probed redshifts. In 
addition, the model predicts an excess of red galaxies at low red- 
shift. It could be argued that part of the discrepancies between 
the SAM and observed colours could be related to uncertainties 
in the modelling of dust extinction (e.g. Kitzbichl er & White| 
|2007[|Fontanot et al.|2009b) . The dashed line in Fig.|7jshows the 
rest-frame B—I colour distribution in the SAM without including 
dust extinction in the model galaxies. In that case, when compar- 
ing to VVDS-Deep colour distributions, one finds that while in 
the highest-redshift intervals the predicted and observed colour 
distributions are quite similar, the lack of blue galaxies in the 
model is remarkable at z < 1.1. Franzetti et al. (2007) show that 
the fixed apparent magnitude selection of the VVDS-Deep sam- 
ple can in principle introduce a mild bias in the intrinsic colours, 
partially displacing rest-frame U — V colours towards the blue at 
z > 1.2. However, they find that this effect is marginal and can- 
not be invoked to explain the discrepancies found in the SAM at 
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Fig. 7. Comparison of the rest-frame B - I distribution in six 
magnitude-limited samples from z - 0.2 to z = 2.1. In each 
panel, the histogram corresponds to the VVDS-Deep measure- 
ment while the curve and associated shaded area correspond to 
the mean and the l<x field-to-field dispersion among the Omocks 
samples. The dashed curves are the mean distribution in the 
Omocks without including dust extinction in the model, while 
the dotted ones are the mean predictions of the Cmocks. 



all probed redshift. Intrinsically, the SAM may not form enough 
very blue galaxies. 

Because of the discrepancies between the predicted and ob- 
served colour distributions, it is difficult and possibly meaning- 
less to define "blue" and "red" galaxies using the same colour 
cut. We therefore opted to use a different colour cut for the SAM 
and the VVDS-Deep sample, with the aim of separating blue 
and red populations on the basis of the colou r bimodality. In the 
VVDS-Deep we use a cut at (B-I) cut = 0.95.|Zuccaetal.|(|2006i 



show that galaxies selected above and below this value largely 
overlap with those classified as early- and late-type galaxies us- 
ing a more refined method based on spectral energy distribution 
fitting. We adopt a larger value of (B - /)"" = 1.3 to separate 
red and blue populations in the SAM. We compare in Fig.|8]the 
total number and the fraction of red and blue galaxies at different 
redshifts. 

We find a large difference in the number density of blue 
galaxies as predicted by the SAM and observed in the VVDS- 
Deep. While the trends with redshift are rather similar, the SAM 
predicts 50-80% more blue galaxies than observed over the 
whole redshift interval 0.2 < z < 1.6 in the VVDS-Deep. In 
the observations we find the presence of an already significant 



number of red galaxies at z — 1.5. This number increases until 
z — 0.8 and then slightly decreases with cosmic time. In contrast, 
the SAM predicts a monotonic increase with time of the number 
of red galaxies. Above z — 0.8, the total number of red galax- 
ies is smaller in the model than in the observations but the trend 
reverses at later epochs, where the total number of red VVDS- 
Deep galaxies starts to decline and that of red model galaxies 
continues to rise slowly. These same problems are evident when 
looking at the fractions of the two populations (bottom panel 
in Fig. |8j. The SAM predicts a monotonic decrease (increase) 
of the fraction of blue (red) galaxies with time at variance with 
VVDS-Deep sample, in which this trend reverses at z < 0.8. 
Note that we discuss the variation with cosmic time of the num- 
ber and fraction of red and blue galaxies in terms of apparent 
increase or decrease in a given apparent magnitude range, which 
does not necessarily imply real increase or decrease in number 
density at these redshifts. 

At redshifts higher than z — 0.8, the number of red galaxies 
is slightly lower in the SAM with respect to observations, while 
blue galaxies appear to be significantly more abundant at all cos- 
mic epochs. This indicates that the overabundance of galaxies 
previously seen in the redshift distribution below z — 0.8 - 1 in 
the SAM (Fig. |5), is due to the presence of a larger number of 
both blue (true at all redshifts for this colour) and red galaxies. 

The difference in the variation with cosmic time of the frac- 
tion of red and blue galaxies and of the shape of the rest-frame 
B — I colour distribution in the SAM, suggests that these popu- 
lations may have different histories of formation and evolution 
than in the VVDS-Deep. In particular, SAM red galaxies may 
start to form at later epochs and be forming continuously and 
more efficiently up to present day, in contrast with what appears 
to be happen to VVDS-Deep galaxies. 

4. Galaxy clustering 

The clustering properties of galaxies provide strong constraints 
on galaxy formation models as they encode important informa- 
tion on how galaxies populate dark matter haloes. We compare 
in this section the galaxy clustering as inferred from the two- 
point correlation function in the SAM with VVDS-Deep mea- 
surements, for the global population of 17.5 < / < 24 galaxies. 

4.1. Two-point correlation function estimation 

We estimate the real-space galaxy clustering using the standard 
projected two-point correlation function, w p (r p ), that corrects for 
redshift-space distortions due to galaxy peculiar motions. This 
is obtained by splitting the galaxy separation vector into two 
components, r p and n, perpendicular and parallel to the line of 
sight respectively (Peebles 1980; Fisher et al. 1994| l, and project- 
ing the two-dimensional two-point correlation function %(r p ,7i) 
along the line-of-sight: 



w p (r p ) = 2 \ %(r p ,n)dn 
Jo 



(1) 



We use the standard [Landy & Szalay ( 1993[ ) estimator to com- 
pute t;(r p ,n). In practice, to obtain w p (r p ), we integrate %{r p ,n) 
up to n max = 20 hT x Mpc. We adopt this value because we find 
that, given the volume of the survey, this value is large enough 
so as to minimise the noise introduced at large n by the uncor- 
rected pairs in the data ( [Polio et al.|2005 1. Errors in the VVDS- 
Deep are estimated through the blockwise bootstrap resampling 
technique (e.g. |Porciani & G iavalisco 2002 ), which allows us 
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Fig. 8. Total number and fraction of red (filled circles) and blue 
(filled squares) galaxies as a function of redshift in the VVDS- 
Deep and in the SAM. The symbols with error bars are the 
VVDS-Deep measurements while solid curves and associated 
shaded areas correspond to the predicted mean fractions and lcr 
field-to-field dispersions among the Omocks. The dotted curves 
are the mean predictions of the Cmocks. In both upper and lower 
panel, higher curves and symbols correspond to blue galaxies. 



to account for sample variance in the field and provides fair er- 
ror estimates, very similar to those obtained using the Jackknife 
resampling technique (Norberg et al.|2009| l. Since the transverse 
dimension of the survey is small, to generate the different resam- 
plings we divide each sample in slices along the radial direction 
(e.g. de la Torre et al. 2009 Meneuxet al.|2Q09]>. For each of our 



samples, we use 300 resamplings by bootstrapping 6 slices of 
equal volume. We estimated the errors in the SAM by comput- 
ing the field-to-field variance among the 100 mock samples. We 
explicitly verified that the bootstrap error estimates agree with 
the ensemble errors obtained from the mock samples. 

The VVDS-Deep sample has a complex angular sampling 
as shown in Fig. [2] To measure the projected correlation func- 
tion we follow the method originally introduced by |Pollo et al. 



2005 ) as improved by |de la Torre et al.| ( |2009| l. This method 
allows us in particular to correct w p (r p ) measurements for the 
inhomogeneous sampling on the sky and the incompleteness on 
small angular scales due to slit mask design. The angular inho- 
mogeneous sampling, i.e. the fact that the sampling rate varies 
with the angular position, is accounted for by reproducing in the 
random sample, used to estimate %(r p ,n), the same variations of 
sampling. This entails accounting for the precise shape of the 
VIMOS field-of-view and the coordinates of the observed point- 
ings. A similar technique was applied in the previous analysis 



of the VVDS-Deep ((Folio et al.||2006 


Meneux et al.||2006 de 


la Torre et al.|2007| Meneux et al.|2008 
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used here includes a more accurate pair weighting scheme to 



Fig. 9. Mean projected correlation function at 0.7 < z < 1.1 for 
red (top curves), all (middle curves), and blue (bottom curves) 
galaxies in the Omocks (dashed curves) and Cmocks (solid 
curves) samples. The error bars correspond to the lcr field-to- 
field dispersion among the mock samples. Omocks points have 
been slightly displaced along r p -axis to improve the clarity of 
the figure. 



correct for missed angular pairs. Indeed, we correct for the in- 
completeness on small angular scales by weighting each galaxy- 
galaxy pair by the ratio of the number of pairs in the spectro- 
scopic sample to that in the parent photometric catalogue (free 
from angular incompleteness), as a function of the angular sep- 
aration. We refer the reader to |Pollo et al.| ( |2005] l and |de la Torre | 
et al. (2009) for the full description of the method used to ac- 
count for the survey selection function and observational biases 
in the measurement of w p (r p ). 

In Fig. [9] we present the estimated galaxy projected corre- 
lation function at 0.7 < z < 1.1 in the mock samples when 
including (Omocks) or not (Cmocks) the detailed VVDS-Deep 
selection function for all, red, and blue galaxies. We find that the 
method does not introduces any systematic errors on the ampli- 
tude and shape of w p (r p ) up to r p = 15 IT 1 Mpc for the blue 
galaxy population. For more strongly clustered population as in- 
carnated by all and red galaxies, our method tends to underesti- 
mate the amplitude of w p {r p ) on scales smaller than ~ 2 - 3h _1 
Mpc by at maximum 15%. In fact, although the method is quite 
robust, we cannot entirely recover the small-scale clustering in- 
formation as we only sample, on average, 27% of the galaxies in 
spectroscopy. However, these systematic errors are smaller than 
the statistical lcr errors of the mock samples. In the following, 
we will then use the Omocks and the VVDS-Deep secure red- 
shift sample to perform clustering comparisons. 

While in general the shape of the observed correlation func- 
tion deviates from a pure power-law (e.g. Zehavi et al.| [2004), 
this simple parametrisation allows one to quantify the clustering 
properties of galaxy samples and to easily compare them. From 
the real-space two-point correlation function £(r), the correla- 
tion length ro, which characterises the clustering strength, and 
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the slope y are obtained by fitting %(r) to a power-law such as 
ij(r) = (r/roY 7 . In the case of the projected (real-space) correla- 
tion function w p (r p ), the power-law form transforms to ( |Peebles| 
[19801 , 



w p (r p ) = r p 



'0 



EfiMg) 



(2) 



where F is the Euler Gamma function. In the present analysis we 
fit the w p (r p ) measurements on the range 0.1 fr 1 Mpc < r p < 
10 fr 1 Mpc using the generalised x 1 method. We use the full 
covariance matrix estimated from the measurements to account 
for the correlations between the different r p bins in w p (r p ) (e.g. 
|Pollo et al.|2005) . 

4.2. Clustering of the global population 

The clustering evolution of the global population of galaxies in 

k_ ((2005b) 



the VVDS-Deep has been measured by Le Fevre et al. 



in the six redshift intervals from z = 0.2 to z = 2.1 previously 
defined. We update their measurements using the more accurate 
method described in the previous section. The w p (r p ) measure- 
ments are shown in Fig. [TO] along with the mean w p (r p ) among 
the SAM Omocks samples. We fit the w p (r p ) with power laws 
and provide the best-fitted parameters ro and y in Tab. [JJ and 
Tab.[2]for both the VVDS-Deep and the SAM. We compare the 
measured correlation lengths in Fig. [TTJ This figure shows that 
our measurements agree very well within the uncertainties with 
those previously obtained by |Le Fevre et al.| ( |2005b[ ) for the same 
field, as well as with results from the DEEP2 survey (Coil et al. 
[2004) 1. 

We recall that for all clustering comparisons, we use model 
predictions computed from the Omocks using the same method 
adopted for the VVDS-Deep data. The Fig. 10 shows that, when 



comparing Omocks (solid curves) and Cmocks (dashed curves) 
measurements, one find a non-negligible bias in the estimation 
of w p (r p ) on sm all scales at z < 0.7. As already pointed out by 



Le Fevre et al.|([2005b|l, here the incompleteness effects are en- 



hanced by the small volume and by the small number of objects 
in the observed samples at these redshifts. This emphasises the 
importance (and the necessity) to include detailed observational 
selection functions and biases in the mock samples in order to 
carry out a fair comparison between observational measurements 
and model predictions. 

At z < 1.1, the SAM predicts on average a higher clustering 
amplitude than measured in the VVDS-Deep, while the slopes 
of the predicted and observed correlation functions are similar. 
In addition, while the correlation length in the VVDS-Deep in- 
creases slightly with increasing redshift, the model predicts a 
roughly constant correlation length over all the redshift range 
probed by the observations. Some evolution of the overall clus- 
tering is expected, because by selecting galaxies at increasing 
redshift in a magnitude-limited sample, we probe intrinsically 
more luminous galaxies on average. Observations show that 
brighter galaxies are more strongly clustered than their fainter 
counterparts (e.g. Norbe rg" st al.|200Tj|Zehavi et al.|2005| Polio 
|et al.||2006) |Coil et al.||2006| l. We~explicitly show in the right 
panels of Fig.|7jand in Tab.[TJand Tab.[2]that we indeed select in- 
trinsically brighter galaxies with increasing redshift, both in the 
VVDS-Deep and the SAM. The relatively constant amplitude of 
the predicted correlation function suggests the absence (or a very 
weak) luminosity dependence of galaxy clustering in the model 
(see also |Li et al.|2007||Kim et aL|2009al l. We study this aspect 
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Fig. 11. Comparison of the correlation lengths as a function of 
redshift. The filled triangles with error bars are VVDS-Deep 
measurements while the dashed curve and associated shaded 
area correspond to the correlation length and error from the 
Omocks. We also report in this figure the previous VVDS-Deep 
measurements by Le Fevre et al. (2005b]) (crosses) and the 
one obtained by Coil et al. ( 2004[ ) in the DEEP2 survey (open 
square). 



in more detail in paper II, where we measure the clustering of 
galaxies with different luminosities in the model and compare 
them with VVDS-Deep observations. 

At z > 1-1, Fig. [TO] shows that, although the overall ampli- 
tude of the predicted galaxy correlation function is similar to that 
measured, the model predicts a shallower correlation function 
than observed. The difference in slope and in shape of the corre- 
lation function can be interpreted within the framework of Halo 
Occupation Distribution models (e.g. |Cooray & Sheth 2002]). In 
this framework, the galaxy correlation function is the sum of 
two contributions, one dominating the smaller scales that char- 
acterises the clustering of galaxies residing in the same halo (the 
1-halo term), and a large-scale contribution, which characterises 
the clustering of galaxies belonging to different haloes (the 2- 
halo term). At z > 1.1, the "bump" of the correlation function 
observed on scales smaller than or of the order of the typical halo 
radius (1 - 21T 1 Mpc) in the VVDS-Deep, suggests that these 
galaxies are on average hosted by relatively massive haloes, with 
relatively large virial radii (e.g. Abbas et al. 2010). The model 
predicts instead a rather weak 1-halo term, which implies the 
presence of relatively few satellite galaxies. Satellite galaxies are 
defined as the galaxies residing within the virial radii of haloes 
that are not associated with their centres. The amplitude of the 
1-halo term is directly linked to the amount of central-satellite 
and satellite-satellite pairs, and in turn to the abundance of satel- 
lite galaxies in their host haloes (e.g 



JBenson etal . 2000; B erlind 
|et al.|[2003] |Kravtsov et al.|[2004| ). We have previously shown 
that at z > 1-1 the model galaxy population is dominated by 
blue galaxies. Therefore, the small-scale shape of the correla- 
tion function predicted by the SAM at these redshift, suggests 
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Fig. 10. Comparison of the projected correlation functions as a function of redshift. In each panel, the filled triangles with error 
bars correspond to VVDS-Deep measurements while the solid curve and associated shaded area correspond to the mean and the lcr 
field-to-field dispersion among the Omocks samples. The dashed curves are the mean Cmocks predictions, while the dotted ones are 
the mean w p (r p ) obtained by rescaling the predicted correlations function in the Omocks to <x 8 = 0.81 as described in Sec. 4.2 



that model galaxies are more likely to be blue central galaxies in 
low-mass haloes rather than blue or red satellites in more mas- 
sive haloes, in contrast with what appears to be the case in the 
VVDS-Deep data. 

It is worth noting that the Millennium Run simulation was 
carried out adopting a WMAP 1 cosmology, with a normalisa- 
tion of the power spectrum of cr 8 = 0.9. The most recent and ac- 
curate measurements favour instead a lower value of <x 8 =s 0.8 1 
(Komatsu et al.|[2009| |20 1 0)> . The value of <x 8 has a non negli- 



gible influence on the amplitude of the dark matter correlation 
function, and in turn to that of galaxies. To quantify the effect on 
SAM clustering predictions, we convert the correlation functions 
in the SAM to those expected assuming the more recent value of 
erg. To do so, we multiply the SAM w p (r p ) by the ratio of the 
non-linear projected correlation function of mass for cr 8 = 0.81 
to that for <x 8 = 0.9. We keep all the other cosmological parame- 
ters identical. We use the |Smith et al.| p003) analytical prescrip- 
tion for the non-linear mass power spectrum, which we Fourier 
transform to obtain the correlation function. With this procedure, 
we assume a purely gravitational clustering evolution and do not 
account for a possible dependence of galaxy bias on <t 8 . This is, 
however, found to be weak in simulations (Wang et al. 2008}. By 
rescaling the SAM projected correlation functions to <x 8 = 0.81, 
we lower the amplitude of the predicted w p (r p ) on all scales and 



vations, there are still some differences, in particular regarding 
the shape of the correlation functions on small scales. 



4.3. Colour-dependent galaxy clustering 

When selecting galaxies at 17.5 < / < 24 in the redshift interval 
0.2 < z < 2, we probe an evolving mix of galaxy luminosities 
and colours (e.g. Le Fevre et al.|2005b Ilbert et al.|20 05 ; Zucca 
|et al. ||2006t |Frarrzetti et al.||2007| l. As a consequence, the clus- 
tering of these galaxies should reflect the average clustering of 
the different galaxy sub-populations at the different redshifts. In 
order to better understand the discrepancies found between the 
global clustering measured from VVDS-Deep data and that pre- 
dicted by the model, we compare the clustering of two sub-sets 
of galaxies, red and blue ones, classified according to their bi- 
modal colour distribution. 

To keep a significant number of galaxies in each redshift in- 
terval, we used only 3 intervals between z = 0.2 and z = 2.1 
in this part of the analysis. We do not consider any sample of 
red galaxies above z = 1.1, as the small number of objects pre- 
vents us from obtaining a robust w p (r p ) measurement at these 
redshifts. Fig. [12] and Fig. 13 show the projected correlation 



obtain the dotted curves in Fig. 10 and Fig. 12 While this tends 
to improve the agreement between model predictions and obser- 



functions and the correlation lengths measured in the SAM and 
in the VVDS-Deep data for red and blue galaxies. The mea- 
surements from VVDS-Deep are consistent with those obtained 
from the same sample by M eneux et al.| ( |2006| >, as well as with 



12 



S. de la Torre et al.: Comparison of VVDS with Munich semi-analytical model high-redshift galaxy properties 



,_, 1000 = 



o 



o 

(X 



100 - 




1 



„ 1000 = 



ioo -i 



10 r 



n — i — i i i i 1 1 1 1 — i — i i i 1 1 1 1 

0.2<z<0.7 



= 1000 



- 100 




i i i 1 1 n il — i i i 1 1 m i 

0.7<z<l.l 



I I I I I III 1 I I I I III 



l.Kz<2.1 



■f-4 



"r-f 



_i i i 



10 



1 

1000 



- 100 



10 



1 



= 1000 



- 100 



= 10 



0.1 



10 



r p [h _1 Mpc] 



s 

J3 3 




SAM Omocks 

• ■ red/blue VVDS (this paper) 
♦ (B-1) >1 3 VVDS (this paper) 

* A red/blue VVDS (Meneux et al. 2006) 
O □ red/blue DEEP2 (Coil et al. 2004) 



0.2 



0.4 



0.6 



0.8 1 
redshift z 



1.2 



1.4 



1.6 



Fig. 13. Comparison of the correlation lengths of red and blue 
galaxies as a function of redshift. The filled symbols with er- 
ror bars are VVDS-Deep measurements while the dashed curves 
and associated shaded areas correspond to the correlation lengths 
and errors in the Omocks. We also report previous VVDS-Deep 
measurements by Meneux et al. (2006 ) as well as those of Coil 
|et al.| ( p004) l from DEEP2 survey, both using different colour cri- 
teria. We refer the reader to the inset for the detail of the plotted 
symbols. 



Fig. 12. Colour-dependent projected correlation function in three 
redshift intervals from z = 0.2 to z = 2.1, both observed in the 
VVDS-Deep and predicted by the SAM. In each panel, the filled 
circles (red galaxies) and filled squares (blue galaxies) corre- 
spond to VVDS-Deep measurements, while the dashed curves 
and associated shaded areas correspond to the mean and the la- 
dispersion among the Omocks. The dotted curves are the mean 
w p (r p ) obtained by rescaling the predicted correlations function 
to cr 8 = 0.81 as described in Sec. 4.2 In all panels, both for 
VVDS-Deep and SAM galaxies, the w p (r p ) with higher global 
amplitude corresponds to that of red galaxies. 



those obtained from the DEEP2 sample by Coil et al. (2004). 
We note that these two studies adopt different colour criteria to 
define blue and red galaxies and that the two surveys have dif- 
ferent observational strategies. In particular, the DEEP2 survey 
selects galaxies brighter than those in the VVDS-Deep, which 
explains the slightly larger observed correlation lengths in the 
DEEP2 survey (see Le Fevre et al.|2005b for a detailed discus- 
sion). 

We find that the correlation functions of blue galaxies in the 
model are in quite good agreement with VVDS-Deep measure- 
ments. In contrast, red model galaxies show a much stronger 
clustering on all scales. This suggests that the stronger cluster- 
ing predicted by the SAM for the entire sample is due to the very 
strong clustering of red galaxies. As for the redshift trend, both 
VVDS-Deep observations and SAM predictions show a rather 
similar evolution over the redshift interval probed. The correla- 
tion length of red SAM galaxies is, however, much larger than 



observed and in fact similar to that observed for extremely red 
objects in the real Universe (e.g., Daddi et al. 2003]) . 

As explained earlier, we have used different rest-frame B — I 
colour cuts to define red and blue galaxies in the model and in 
the observations. One may argue that the higher clustering am- 
plitude observed for SAM red galaxies could be due to the redder 
colour cut applied to select the two populations, redder galaxies 
being expected to be more strongly clustered. To test this possi- 
bility we measure the clustering strength of red galaxies in the 
VVDS-Deep using the same colour cut used in the SAM, i.e. 
(B - I) cl " = 1.3. In this way, we isolate in the VVDS-Deep sam- 
ple red galaxies which have the same rest-frame B - I colour 
distribution than red model galaxies. As shown in Fig. 
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find that these galaxies in the VVDS-Deep show a higher clus- 
tering amplitude than those selected with B-I > 0.95. However, 
the SAM clustering remains significantly stronger, demonstrat- 
ing that red model galaxies are intrinsically more clustered than 
observed in the VVDS-Deep. 

When studying the shape of the projected correlation func- 
tions in more detail, one finds that blue SAM galaxies are char- 
acterised by a shallower correlation function than VVDS-Deep 
galaxies, in particular on small scales. Within the HOD frame- 
work, this implies a weaker 1-halo term that can be interpreted 
as a lack of blue satellite galaxies in the SAM. In contrast, red 
model galaxies exhibit a correlation function which is signifi- 
cantly steeper and higher than observed. Here, the very promi- 
nent 1-halo term may be due to an overabundance of red satellite 
galaxies. Similarly, Coil et al. (2008 ) find an absence of "Finger 
of God" (FoG, Jackson 1972) in the correlation function of blue 
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model galaxies at z — 1 at variance with red model galaxies, 
which have a very strong FoG. The FoG effect is associated with 
the infall of satellite galaxies inside haloes and its strength is 
related to the abundance of satellite galaxies (e.g. |Slosar et al.| 
2006). These results suggest that in the real Universe, (at least 
part of) the red satellites likely evolve less rapidly than in the 
model, and remain in the blue tail of the colour distribution for 
a longer time scale. This could adjust the different small-scale 
clustering behaviours of blue and red SAM and VVDS-Deep 
galaxies, but would not affect significantly the amplitude of the 
correlation functions. 

To better see the impact of an overabundance of satellites on 
the clustering of galaxies, we randomly remove from the SAM 
mock samples 80% of the red satellites. The resulting galaxies 
correlation functions are shown in Fig. 14 This figure shows 



that, by excluding most of SAM red satellites, the amplitude of 
the correlation function of red galaxies is dramatically reduced, 
particularly on small scales. Model predictions obtained exclud- 
ing 80 per cent of the red satellites are in quite good agreement 
with observational measurements but at 0.2 < z < 0.7, where 
there is still a significant difference between the amplitudes of 
the predicted and measured correlation functions. Similar con- 
clusions have been reached while comparing the model to local 
measurements, e.g. Li et al. ( 2007) 1 found that the match of the 
observed clustering in the local Universe to the previous version 
of the Munich model ( |Croton e"t"a l. 2006) can be improved by 
removing 30% of satellite galaxies in the model. 

5. Summary and discussion 

We have compared some of the basic high-redshift galaxy prop- 
erties as measured in the VIMOS-VLT Deep Survey, to predic- 
tions from the Munich semi-analytical model. For this purpose, 
we have constructed 100 mock samples that accurately mimic 
the VVDS-Deep observational strategy. We have compared the 
magnitude counts, redshift distribution, colour bimodality, and 
galaxy clustering for galaxies with 17.5 < I < 24, probing a 
broad range of cosmic epochs from z=2 to z — 0.2. We have 
demonstrated that, in order to carry out a fair comparison be- 
tween model predictions and data, it is important to build "ob- 
served" mock samples that accurately reproduce the detailed se- 
lection function and biases of the observations. 

We find that the Munich semi-analytical model reproduces 
reasonably well: 

- the magnitude counts in the u*, g', r', i', I and rest-frame B 
bands, 

- the shape of the redshift distribution at z < 1.8 for Iab < 24 
galaxies, given the relatively large sample variance predicted 
by the model in the VVDS-Deep volume, 

- the global galaxy clustering at z > 0.8, 

but fails to reproduce: 

- the magnitude counts in the z' and rest-frame / bands, 

- the shape of the redshift distribution at z > 2 for Iab < 24 
galaxies, 

- the rest-frame B - 1 colour distribution and its evolution with 
cosmic time since z — 2, 

- the clustering strength of red galaxies, 

- the detailed small-scale clustering of both red and blue 
galaxies. 

It is important to notice that for some of the predicted galaxy 
properties, there is a significant variance among different mock 
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Fig. 14. Red and blue galaxy projected correlation functions in 
two redshift intervals from z = 0.2 to z = 1.1, both observed in 
the VVDS-Deep and predicted by the SAM (mean over Omocks 
samples). In each panel, the solid curves correspond to SAM 
mean predictions while the dashed ones to the resulting mean 
predictions while keeping only 20% of red satellite galaxies. 
The filled circles (red galaxies) and filled squares (blue galaxies) 
correspond to VVDS-Deep measurements. Both for the VVDS- 
Deep and the SAM, the w p (r p ) with higher amplitude corre- 
sponds to that of red galaxies. The SAM correlation functions 
are rescaled to <x 8 = 0.81 as described in Sec.|4.2| 



samples. For most of the observational measurements discussed 
in this study, we find that there are a few mock samples that are 
in good agreement with the data. On average, however, models 
deviate from observational measurements. In particular, for the 
colour distribution and the clustering of red galaxies, all mock 
samples differ from the VVDS-Deep measurements, and differ- 
ences are larger than 3cr. None of the mock samples is able to 
reproduce all the VVDS-Deep measurements presented in this 
analysis, suggesting that the model failures highlighted above 
are not simply due to sample variance. 

The discrepancies found between model predictions and 
VVDS-Deep observations extend to higher redshifts some of 
the model problems that have been previously emphasised from 
data-model comparisons in the local Universe. Although the 
blue population dominates in number density at all redshifts, the 
SAM tends to produce too many relatively bright red galaxies. 
As a consequence, the rest-frame /-band distribution is skewed 
towards bright magnitudes and the rest-frame B - I colour dis- 
tribution towards the red. This excess of red galaxies is dom- 
inated by satellites, giving rise to a prominent 1-halo term in 
the correlation function of red model galaxies. In addition, the 
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SAM underpredicts the fraction of blue satellite with respect to 
blue central galaxies as seen in the small-scale clustering of blue 
galaxies. It is important to mention that the excess of red satellite 
galaxies is not specific to the Munich semi-analytical model but 



is present in most of published semi-analytical models (Liu et al. 
[2OT0l l. 

The excess of red and deficit of blue satellite galaxies 
in semi-analytical models are likely due to an over-efficient 
quenching of satellites, that transforms too many blue galaxies 
to red ones over a short time scale. In the models, star-forming 
blue galaxies become passive as a consequence of the infall of 
a galaxy onto a larger halo, or because of AGN feedback that 
suppresses star formation in massive central galaxies. An over- 
quenching of satellite galaxies can be produced by a too efficient 
strangulation, that instantaneously shuts off the star formation 
when a gal axy enters in a halo ([ Weinman n et al.||2006| |Font| 
et al.||2008| |Kang & van den Bosch||2008| |Kimm et al.||2009| 
Fontanot et al. 2009a). It is interesting to note that the trunca- 



tion of gas accretion in satellite galaxies, and to a large extend 
of star formation, is also found to be less abrupt in smoothed 
particle hydrodynamics simulations (Cattaneo et al. 2007 ; [Saro| 
et al.|[20T0] >. This would help to explain the difference in the 



rest-frame colour distribution between models and data, but can- 
not explain the very strong intrinsic clustering of red galaxies. 



As recently pointed out by Kim et al. (2009a), Wetzel & White 
(2010), and |Liu et al] ( |2010[ ), fheprobl em might lie in a poor 



treatment of satellite mergers and disruption. In fact, most of 
current semi-analytical models, including that used in this study, 
do not account for tidal stripping of satellite galaxies (but see 
[Benson et al.|2 002; Monaco et al. 2007). This can influence sig- 
nificantly the predicted clustering signal, and has been shown to 
affect also the predicted galaxy intrinsic colour distribution and 
stellar mass function ( |Yang et al.|20 09 ). We study these aspects 
in more details in paper II where we specifically relate the pre- 
dicted galaxy clustering as a function of luminosity and colour 
to the halo occupation predicted by the model. 
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